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Abstract. We study the isochronicity of centers at O G M for systems 

±^ -y + A{x,y), y = x + B{x,y), 

where A, B £ M.[x, y], which can be reduced to the Lienard type equation. When deg{A) < 4 
^ , and deg{B) < 4, using the so-called C-algorithm we found 36 new multiparameter families of 

QQ ' isochronous centers. For a large class of isochronous centers we provide an explicit general formula 

^H for linearization. This paper is a direct continuation of 1 3 1 but can be read independently. 

o' 
in 

i/^ . 1. Introduction 

^— V ' Let us consider the system of real differential equations of the form 

dx A/ \ dy n / \ 

— == X = -y + A(x,y), — ^ y ^ x + B(x,y), (1.1) 

where (x, y) belongs to an open connected subset [/ C K^ containing the origin O = (0, 0), with 
A, B e C^{U, M.) such that A and B as well as their first derivatives vanish at O. An isolated singular 
C^ I point p £ U of system dl.ll ) is a center if there exists a punctured neighborhood V £1 U of p such 

that every orbit of dl.lb lying in y is a closed orbit surrounding p. A center p is isochronous if the 
period is constant for all closed orbits in some neighborhood of p. 

The simplest example is the linear isochronous center at the origin O given by the system 

X = -y, y^x. (1.2) 

The problem of characterization of couples [A^ B) such that O is an isochronous center (even a 
center) for the system ( 11.11 ) is largely open. 

The well known Poincare Theorem asserts that when A and B are real analytic, a center of ( 11.11 ) 
at the origin O is isochronous if and only if in some real analytic coordinate system it takes the form 
of the linear center ( 11.2b (see for example [1], Th.13.1, and |20|, Th.4.2.1). 

An overview 1,6] presents the basic results concerning the problem of the isochronicity, see also 
ifTl [TOl fm [121 I20II ■ As this paper is a direct continuation of |J3l, we refer the reader to it for general 
introduction to the subject. Here we will recall only the strictly necessary facts. 

In some circumstances system dl.lb can be reduced to the Lienard type equation 

x + f{x)x^ +g{x) = Q (1.3) 
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with f,gE C^{J,M.), where J is some neighborhood of G K and g{0) — 0. In this case, 
system ( II. Il l is called reducible. To equation ( 11.3b is associated the equivalent, two dimensional 
(planar), Lienard type system 

x = y 

y = -9{x) - f{x)y' ^ 

For reducible systems considered in this paper, the nature (center and isochronicity) of the singular 
point O for both systems ( 11.11 ) and ( 11.4b is the same. 

Let us return now to the Lienard type equation ( 11.3b . Let us define the following functions 



2 r ■ (1-4) 



F{x):= fis)ds, 0(x) := / e^'-^'ds. (1.5) 

^0 Jo 

The first integral of the system (11.4b is given by the formula (f2Tl, Th.l) 

I{x, x) = ^(ie^("))' + r .g(s)e2-P'(^)ds. (1.6) 

2 Jo 

When xg{x) > for x ^ 0, define the function X by 

\a^f ^ f 9{s)e^^^'^ds (1.7) 

and x^{x) > for a; 7^ 0. 

Theorem 1.1 (EH, Theorem 2). Let f, g eC^ (J, M). Ifxg{x) > Ofor x ^ 0, f/ien f/ie system ( fTH l 
/jfli fl center at the origin O. When f and g are real analytic, this condition is also necessary. 

Theorem 1.2 (| 12], Theorem 2.1). Let f and g be real analytic functions defined in a neighborhood 

JofO € M., and let xg{x) > Ofor x ^ 0. Then system ( 11.4b has an isochronous center at O if and 
only if there exists an odd function h G C^{J, M.) which satisfies the following conditions 



^i^) /^\.^F{x 



1 + hia^)) 



gix)e'^^>, (1.8) 



the function (t>{x) satisfying 



(l){x) ^i{x)+ / h{t)dt, (1.9) 

JQ 

and £,{x)(f>{x) > Oforx 7^ 0. 

When those equivalent conditions are satisfied, then the function h is analytic in the neighbor- 
hood J. It is called the Urabe function of system (11.4b . 

Taking into account ( 11.71 ). it is easy to see that ( 11.8b and ( 11.91 ) are equivalent. 

Corollary 1.3 (| 12], Corollary 2.4). Let f and g be real analytic functions defined in a neighbor- 
hood 0/0 € R, and xg{x) > Ofor x ^ 0. The origin O is an isochronous center of system ( 11.4b with 
Urabe function h = if and only if 

g\x)+gix)f{x)^l (1.10) 

for X in a neighborhood of 0. 

In the sequel we shall call the Urabe function of the isochronous center of reducible system 
dl.lb the Urabe function of the corresponding Lienard type equation. 

In IT2I the third author described how to use Theorem ll.2l to build an algorithm (C-Algorithm, 
see Appendix of ID for more details) to look for isochronous centers at the origin for reducible 
system dl.lb . and applied it to the case where A and B are polynomials of degree 3. This work was 
continued in I.13J and in L3J. 
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In this paper we apply the so called Rational C-Algorithm introduced in fT\ which is more 
powerful then the algorithm used in fSl . 

The aim of the present paper is to extend these studies to the following real multiparameter 
family of polynomial system of differential equations : 



i = -y + ai,ixy + a2,QX^ + a2,ix^y + a^^x^ + a^^ix^y + ai,Qx'^ \ 

y = X + bo.2l/ + bi.ixy + b2,ox^ + bi_2xy^ + b2.ix^y + bsxtx^ + 62,22;^?/^ + bs^ix^y + 64,02;'* J 

(1.11) 
The reported results, which are obtained by Maple computations are reproduced without almost any 
change to avoid misprints. 

The paper is organized as follows. In Section |2] we report the necessary background and de- 
scribe the investigated subfamilies of system (II. lib . In Sections [3] and H] we describe the obtained 
new isochronous centers. In total we provide 36 new families of isochronous centers. Among them 
two Monsters (14.23b and (14.30b of extreme complexity, never encountered before. 

Let us stress that when describing the Urabe functions of the isochronous centers from Section 
3, for the first time we encounter the non-standard examples of it. Indeed, up to now all identified 
Urabe functions was always of the form h(f) — "^ where a, b, c £ M., b > and n a non 

negative integer.(See lfT2llT3l l3l') 

Finally, in Section |5] when Urabe function /i = 0, we describe the explicit general formula 
for linearizing change of coordinates whose existence is insured by the Poincare theorem. We report 
also 5 examples of such linearization. 



2. Preliminaries 

2.1. Choudhury-Guha Reduction 

Let us consider the real polynomial system 



X ^po{x) +pi{x)y 

y ^qo (x) + qi {x) y + q2 (x) y 



2 r ' (2-1) 



wherepo,Pi,go,'7i,<72 e R[a;]. 

We will always assume that O = (0, 0) e M^ is a singular point of ( 12.1b . that is po{0) 
9o(0) = 0. Let us assume also thatpi(O) 7^ 0. 

Let us note that the system (II. lib is a particular case of (12.1b when 

po (x) = a2^ox^ + aa^ox^ + 04,02;'' 

Pi (x) = -1 + fli^ix + 02,12;^ + 0342;^ 

qo (x) = x + b2fiX^ + &3,q2;^ + 64,02;'* ) ■ (2.2) 

qi {x) = 614.T + 62.12;^ + 63, ix^ 

<?2 (2;) = 60,2 + 61,22; + 62,22;^ 
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The following change of coordinates x ~ x, z = po {x) +pi (x) y transforms the system ( 12. Il l 
to the system 



X = z 



92 {^) , P[ {x)\ ^, ^ f{p'Ax))p,{x) ^ ^^ ^^^ _^^, ^^^ _ 2 ^iMznM , , 



pi{x) Pi{x)J \ pi{x) pi{x) 

92 (x) {po (x)) 



Pi (x) 



- 91 (x) Po (x) + pi (x) qo (x) 



(2.3) 



If 



[P'l (x)) Po (x) , q2{x)po{x) 

+ 91 (X) + Po (X) - 2 = 0, (2.4) 



Pi (a:) Pi (a;) 

the system ( 12.1b is of Lienard type (ll.4l l. 



with 



92 (a:) ^ Pi (a;) 
pi (x) pi (x) 



;(a;) = - ^ 'y' + qi (x) po (x) - pi (x) qo (x) 
Pi[x) 



(2.5) 



To the best of our knowledge, the above reduction of the system ( 12.1b to Lienard type system (11.4b 
was proposed for the first time in a preliminary and never published version of ||9l. It is then natural 
to name it Choudhury-Guha reduction. 

In all considered cases (see ( 12.21) ) it is easy to see that for |a;| small enough g{x) — x + x'^ g{x) 
where ^ is a real analytic function. Thus xg{x) > for x ^ 0, \x\ small enough and Theorem ll.2l 
insures that the origin O is a center for the system ( 11.4b . Our aim is to decide when this center is 
isochronous. 

When Po and qi identically vanish, ( 12.4b is satisfied and we recover the standard reduction 
from |[3l (see Case 1 from Sec. 1 of [3]). Many particular cases of system dl.lll ) where studied, 
using this standard reduction. 

1. In [12| : aij = 02,0 = 03^0 = 03,1 = 04,0 = ^2,2 = &4,o = &3,i = h.i = &i,i = 0. That 
means thatpo(a;) = qi{x) = 0, pi{x) = — 1 + 02, ix^ and deg(go(a;)) < 3, dcg(q2(a;)) < 1- 

2. In IIT3I : 02,0 = 03,0 = ^3,1 = 04,0 = ^2,2 = &4,o ~ &3,i = ^2,1 — &i,i = 0. That means that 
Po{x) = qi{x) — and we consider only cubic systems. 

3. In lO three families are studied : 

(a) 02,0 = 0-3,0 = 04,0 = ^3,1 = ^2.1 = ^1.1 = with zero Urabe function. 

(b) ai,i = 63,0 = a2,o = 03.0 = a4,o = &3,i = ^2,1 = h,! = 0. 

(c) ai,i = 02,1 = 02,0 = 03,0 = a4,o = ^3,1 = ^2,1 = ^1,1 = 0. 



2.2. Investigated families 

The exhaustive study of all isochronous center at the origin for the system ( 11.111 ) is hopeless at the 
present. Even for cubic system when all quartic terms vanish this problem is not yet solved. 

Let us note that the condition (12.4b is equivalent to the following system of equations: 
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2 02,0 + ^14 = 0, 

- &2,i + ai,i5i4 + 01,102,0 - 2 60,202,0 - 3 03,0 = 0, 
01,162,1 - 4 Gifi -63,1 + 2 01,103,0 - 2 60,203,0 - 2 61,202,0 + 02,161,1 == 0, 

- 03,102,0 + 3 01,104,0 + 02,162,1 - 2 60,204,0 - 2 62,202,0 + 02,103,0 + 01,163,1+ \ . (2.6) 

- 261,203,0 + 03,161,1 = 0, 

03,162,1 + 02,163,1 + 2 02,104,0 - 2 61,204,0 - 2 62,203,0 = 0, 

03,163,1 - 262,204,0 + 03,104,0 = 0. 

In the present paper we determine all isochronous centers of the system dl.llb in each of the 
following three cases. 

1. When the standard reduction is possible (i.e. pq{x) = and qi{x) — 0, that is 02.0 — 03.0 — 
04. = 63.1 = 62,1 = 61,1 = 0). We provide all candidates for isochronous centers in the cases 
where either ai,i = 0, or 62,0 = —3 60,2- In all the cases but one (a subcase of 62,0 — —3 60,2), 
we prove the isochronicity. The general case is not yet completely explored. 

2. When Choudhury-Guha reduction for the cubic case is possible (i.e. conditions ( 12.61 ) are satis- 
fied and 03,1 — 04,0 — 63,1 = 62,2 = 64,0 = ). For this case we obtain the exhaustive list of 
all isochronous centers at the origin. 

3. When Choudhury-Guha reduction is possible and the Urabe function is null. That means that 
condition ( 12.6b and condition dl.lOb for / and g defined by (12.5b are simultaneously satisfied. In 
this case we provide 25 examples of new isochronous centers and our analysis is not exhaustive. 

Moreover, when the Urabe function /i = 0, we give the explicit formulas for linearizing coor- 
dinates from Poincare Theorem. We report 5 examples where such coordinates are explicitly com- 
puted. 

2.3. Time-reversible systems 

The general notion of time-reversible system of ordinary differential equations goes back to ifTSJI 
where the motivations and general discussion can be found. Here we follow Sec. 1 of 15] (see also 
Sec. 3.5 |20|). 

The planar system ( 11.11 ) of ordinary differential equation is time-reversible if there exists at 
least one straight line passing through the origin which is a symmetry axis of the phase portrait of 
the system under consideration. By appropriate rotation this straight line is mapped on the x-axis 
and the phase portrait of the rotated system is invariant with respect to symmetry {x, y) — > (x, —y) 
if only one change time t into —t. 

Note that a system 

X = P{x,y), y = Q{x,y), 

istime-reversiblesystemwithrespecttox-axisif andonly if P(x, — y) — -~P{x,y) andQ{x,—y) — 
Q{x.y). When P and Q are polynomials, this means that the variable y appears in all monomials of 
P in odd power and in all monomials of Q in even power (0 included). 

Consequently, to decide if a polynomial center for system ( II. lb is time-reversible or not, we 
consider the rotated system in coordinates (xcJ/q), where Xq = xcosa—ysina?iViAya — x sina+ 
y cosa and we examine the parity of the powers of the variable y^ for all angles a. 

This notion plays an essential role in our topics. Indeed, for system dl.lb the origin is either a 
center or a focus. Thus, if such system is time-reversible the focus case is excluded and the origin is 
necessarily a center 

To the best of our knowledge the majority of already known isochronous centers for polynomial 
system dl.lb are time-reversible. For instance, all systems studied in ll22l l5l [121 [131 [3l fSl are time 
reversible. Moreover, among 27 polynomial isochronous centers presented in tables 3 — 29 of Q 
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only 7 are not time-reversible; indeed, those from tables 17 and 23 — 28. In what concerns the cubic 
isochronous centers for system ( II. lb the complete enumeration of those which are time reversible 
was obtained in [81; there are 17 such cases. In |6 1 one find 4 non time-reversible isochronous centers 
(tables 25 — 28) and in the present paper we present three new such cases which are described in 
Theorem 14.11 In [7 1 a complete list of quartic homogeneous time-reversible isochronous centers 
is provided, there are 9 such cases. In the present paper we provide 33 new cases of quartic (non 
homogeneous) isochronous centers. Among them 8 are time reversible and at least 23 are not time 
reversible. 

2.4. Background on Grobner bases 

The use of the Rational C- Algorithm leads to a system of polynomial equations 

/i=0,...,/,„ = (2.7) 

with /,; e M[xi, . . . , Xn]- To solve this system we consider the ideal (/i, . . . , /,„) C R[xi, . . . , a;„]. 
For this aim, we use Grobner Bases computations. In this section, we recall the basic facts about 
Grobner bases, and refer the reader to [14] for details. 

A monomial ordering is a total order on monomials that is compatible with the product and 
such that every nonempty set has a smallest element for the order The leading term of a polynomial 
is the greatest monomial appearing in this polynomial. 

A Grobner basis of an ideal I for a given monomial ordering is a set G of generators of I such 
that the leading terms of G generate the ideal of leading terms of polynomials in I. A polynomial 
is reduced with respect to the Grobner basis G when its leading term is not a multiple of those 
of G. The basis is reduced if each element g G G is reduced with respect to G \ {g}. For a given 
monomial ordering, the reduced Grobner basis of a given set of polynomials exists and is unique, 
and can be computed using one's favorite general computer algebra system, like Maple, Magma or 
Singular. The most efficient Grobner basis algorithm is currently F4 ifTSl . which is implemented in 
the three above cited systems. For our computations, we use the FGb implementation of F4 available 
in Maple L17J. 

The complexity of a Grobner basis computation is well known to be generically exponential 
in the number of variables, and in the worst case doubly exponential in the number of variables. 
Moreover, the choice of the monomial ordering is crucial for time of the computation. 

The grevlex ordering is the most suited ordering for the computation of the (reduced) Grobner 
basis. The monomials are first ordered by degree, and the order between two monomials of the same 
degree Xa — x^^ ■ ■ ■ a;"" and xp — x^ ■ ■ ■ x^" is given by Xa >~ xp when the last nonzero element 
of (ai — /3i, ...,«„ — /3„) is negative. Thus, among the monomials of degree d, the order is 

xi > xf'^X2 >- xf~'^xl y ■■■ y X2'^ X'^'^X3 >- xf~^X2X3 >- x'^'^xj ^ • • • )- xf^. 

However, a Grobner basis for the grevlex ordering is not appropriate for the computation of the 
solutions of the system (I2.7l i. The most suited ordering for this computation is the lexicographical 
ordering (or lex ordering for short). The monomials are ordered by comparing the exponents of the 
variables in lexicographical order. Thus, any monomial containing xi is greater than any monomial 
containing only variables X2, ■ ■ ■ ,Xn. 

Under some hypotheses (radical ideal with a finite number of solutions, and up to a linear 
change of coordinates), the Grobner Basis of an ideal (/i, . . . , /,„) for the lexicographical order 
Xi > . . . > Xn has the shape 

{Xi - gi{Xn),X2 - g2iXn),---,Xn-l - 5„_1 (x„_i), g„ (x„)}, (2.8) 

where the (7i's are univariate polynomials. In this case, the computation of the solutions of the system 
follows easily. In the general case, the shape of the Grobner basis for the lexicographical ordering is 
more complicated, but it is equivalent to several triangular systems for which the computation of the 
solutions are easy. 
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An important point is that a Grobner basis for the lex order is in general hard to compute 
directly. It is much faster to compute first a Grobner basis for the grevlex order, and then to make a 
change of ordering to the lex order. 

The precise ordering we use to compute the Grobner bases of the polynomial systems oc- 
curing in this paper is a weighted order: we fix a weight i + j — I for the variables Uij and bij 
(see (II. ni l), and use the weighted grevlex or lex ordering. For those orderings, the polynomials are 
homogeneous, which simplifies the computation. Indeed, without loss of generality, we can pick 
a variable a.ij and split the computations into two cases Uij — and Uij = 1 (the same concerns 
bij). The entire set of solutions can then be recovered in the standard way. For instance, all solutions 
with ai^i y^ for system ( 11.111 ) are obtained from solutions with aii = 1 by the change of variables 
X = Oil a: and Y = ai i y. This trick reduces by one the number of variables for the Grobner basis 
computation and improves the time of the computations. In what follows, all results are presented 
up to such homogenization. 

Finally, we use repeatedly the Radical Membership Theorem: 

Theorem 2.1 ( 11411 ). Let I = (/i , . . . , fs) be an ideal ofk[xi , . . . , x„], then f belongs to y/ if and 
onlyif{fi,...,fsA-yf) = (1) = k[xi,...,Xn,y]- 



3. The standard reduction 

In this section we are concerned by system dl.lll ) with 02,0 = 03,0 ~ ^4,0 = ^3,1 = ^2,1 ~ &i,i = 
which gives : 

X = -y + ai,ixy + a2,ix^y + a3,ix^y 1 

y = X + b2fix'^ + b^fiX^ + 60^2?/ + bi^2xy'^ + b2,2x^y^ + &4,oa;^ J 

Recall that all cases when the origin O is an isochronous center of the system ( 13.11 ) with zero 
Urabe function are described in |3|. In the following theorem we omit all isochronous centers with 
zero Urabe function, as well as all cubic isochronous centers that where all described in llT3irT2ll . 

For each case we prove the isochronicity by determining explicitly its Urabe function. For 

system (k) we will denote it by h(i.y 



Theorem 3.1. The following particular cases of system ( 13.11 ) have an isochronous center at the 
origin O. 



3x^y± V2x^y 

JT) ^F) ) 1 where h^^ 

V^_2 , __3 , A _2 I r, /;:__2_2 I V^__4| ' — ' 



±e 



y^x±V2x'T^y'+x'+ixy^±2V2xY±^x*\ ^2 + 9^ 



X = —y + x'^y 




3x^ y"^ bx^ 

X — — y + 9 x''y + 6 x y + xy 
y^x 



^4 + 49^2 



(3.2) 
(3.3) 

(3.4) 

(3.5) 
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i = -y - (3 03,1 + -] x'^y + a^^ix^y + xy 
y^x+\ -3 034 + - ) xy- 



X = — y + xy 
y^x 



3 X- 



y"^ +x^ 



> , where h^^ = 



, where h^j^ = 



v/(l-27a3,i)e + 9 



(3.6) 



(3.7) 



2 " 4 

Moreover, all other possible isochronous centers at O for non cubic system ( 13.11 ), where either 
ai^i — I or 62,0 — ~3 60.2, and with non vanishing Urabe functions, belong to the family 



-y 



- - 2 ^2,2 )^^y+ [y^ + ^2,2 ) x^y + xy 



y = X 



3x-^ 



2„,2 



2b2,2xy +b2^2x y 



(3.8) 



3^ y2 

^- + J ■ 8 ---'-"- ■ "-^-" - 16. 

In particular, when 62,2 G {^TS' '^i tki' ^'^^ origin O is an isochronous center with non- 
standard Urabe functions: 



V2^ 2L ^- 



Hb2.,=-^} 



e 



L i±'. 



3 L 



<^{ 



f +* M^ 



where L — LambertW is the Lambert function (see [16]), 

V2 



^{62.2=0} 



-4+;^+2^4+2g3 , 



e e' + 2V4 + 2|2 + 2 



^''2.2 = 1^} 



(2 + ^2)(^y4^2e^+6) 

V2gV2e^+32(C^ + 12) 
2 (^2 + 4) ie + 16) 



Proof. Necessary conditions. This part of the proof is based on the C-Algorithm. Indeed, 19 steps 
are necessary to find the algebraic conditions of isochronicity (see Appendix A of |3|). We did not 
succeed in computing the full grevlex Grobner basis of the corresponding system of polynomial 
equations. We restricted ourselves to the cases aii ~ 0, which gives the cases (I3.2l )-( l3.3b . and 
{ai.i — 1, 62,0 = —360,2} which gives the cases ( I3.4b -( l3.6l l and ( 13.81 ). We also get case ( 13.7b as a 
particular solution. 

Sufficient conditions. For each case we determine its Urabe function. For systems ( I3.2I )-( |3^ 
the procedure in Section 2 of |3| is applied. 

The search of the Urabe function for sytem ( 13.8b is more subtle. Indeed, we verified that for 
all values of parameters the first 20 necessary conditions of isochronicity given by C-algorithm are 
satisfied. This strongly suggests that for all values of parameters the system ( 13.81 ) has an isochronous 
center at the origin O. For this system 

20 - 96 0:62 2 + 64 a;2 62 2 - 12 a; + 3 x'^ 



fix) 



-16 - 6 a;2 - 32 a;2&2 2 + 16 x^b2 2 + a;3 + 16 a; 



and 



g{x) = — - (-16-6a;2-32a;262,2 + 16x^62,2+a;''' + 16x)a;(-16 + 12a:--6x2 



256 
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from formula ( 11.71 ) one obtains 

^^(x) = 2x2 (2 - a;)-2(i6fc2,2+i)-^ (4a;262,2 + 1/4x2 -x + 2y'^^^ . (3.9) 

From formula ( 11.8b one deduces that 



h{^{x)) = - 



(x - 4) (x2 - 2 x + 4) ' 



Now the problem is to find the reciprocal function x = x{^). Unfortunately we succeeded in finding 
it only for 52,2 € {~Tf,^ ^^ iki because in those cases the equation ( |3.9l l takes a sufficiently simple 
form. 

D 

Note that the system dS.BI l was already identified in Theorem 2.2 of Q. 



4. The Choudhury-Guha reduction 

4.1. Cubic isochronous centers 

Choudhury-Guha reduction is more general than the standard one used in preceding papers lfT2l[T3l 
O. Here we provide the complete enumeration of all cubic systems from dl.lll) and we find three 
new cases of isochronous centers at the origin. The system we consider is: 



(4.1) 



(4.2) 



X = -J/ + ai^ixy + a2flx^ + a2,ix^y + 03,02;^ 
y ^ x + b2.ox'^ + bi,ixy + bo^2y^ + b2syx^ + h,2y^x + ba^x^ ^ 
Condition ( I2.4l i is equivalent to the following system of equations: 

2 02,0 + 61,1 = 

01.102.0 - 3o3,o - ^2,1 + 01,161,1 - 260,202,0 = 

01.162.1 - 2 60,203,0 + 2 01,103,0 + 02,161,1 - 2 61,202,0 = 
02,103,0 + 02,162,1 - 261,203,0 = 

Theorem 4.1. Under the assumptions ( I4.2l l the origin O is an isochronous center of System (14.1b 
only in one of the following cases: 

1. The standard reduction is possible, that means 03.0 = 02,0 = 61,1 = 62.1 = and the system 
is one of those from Theorem 3 of II13I . 

2. We are in one of the following cases : 

X ^ -y~2b2fixy + x"^ + 2b2,ox^ 1 

y = x-4: 62,02/^ -2xy + b2fix^ + 4 62,oa;2?; + 2 x^ J 



(4.3) 



X = -y ± 2 \/2xy + x^ ip 2 V2x^ 

y = X ± 8 \/2y2 - 2 xy T 3 v^x^ ^ 12 ^/2x^y + 10 x^ 



i = -y - 2 b-ifixy + x^ + - b2flX^ 

y = X - 62,02/2 - 2 xy + 62, oa;^ + b2flX^y + ( 2 + - 62,0^ ) a;^ 



(4.4) 



(4.5) 
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Proof. The necessary conditions are given by the solutions of the polynomial system of equations 
consisting of equations ( 14.2b (called Ci , . . . , C4) and the 8 equations obtained from the Rational 
C-Algorithm (15 steps), called C5, . . . , Ci2- Let us denote by / the ideal generated by Ci , . . . , C12. 
We exclude the standard reduction by adding to / the variable T and the polynomial 

Ci3 = [Ta^fi - l){Ta2fl - l)iTh,i - l){Tb2,i - 1). 

For 02,0 = 0, a Grobner basis of (Ci, . . . , C12, C13, 02,0) is (1) (i-e. there is no solution), which 
implies that we can take 02,0 — 1- We use the weighted order 61 1 > 62,1 > 63,0 > ^1,2 > 02,1 > 
^0,2 > ai,i > &2,o > 03,0- 

First, a Grobner basis of system (Ci , . . . , Cg) for the weighted lex order contains the polyno- 
mial 

P = (ai.i + 2 5o,2)(a2,i - ai, 103,0 - 03,0)- 
We spUt our problem into two subcases according to this factorization. 

• for ai 1 + 2 5o,2 = 0, we get only one real solution 

X ~ —y + X I 

y = X - 2xy + 2x^ J ' 

which is a particular case of ( 14.51 ) with 62,0 = 0. 

• for 02.1 — 01.103.0 — ai = 0, we eliminate the solutions that are not real by adding to the 
system the polynomials Pi ■ Ti — 1 for each P, in 

{16 alo + 1,4 alo + 9, 4 alo + 1, a|o + 4, a|o + 1, 03,0 + 16, 03,0 + 9-4 62,0 03,0 + 4 6^,0} 
that have no real solution. Then, the solutions of the resulting system are those quoted in the 
theorem. 

Sufficiency. For the cases (14.3b and (14.4b we have g'{x) + f{x)g{x) = 1. Hence by Corol- 
lary [T]3]the origin is an isochronous center. Moreover we easily check that /i|43j(^) = — ^62.oC- n 

4.2. Quartic isochronous centers 

Our first target was to identify all isochronous centers at the origin with zero Urabe function for the 
system dl.UI) 

x = -y + ai_ixy + 02, oa^^ + ai.ix^V + 03,02;^ + aj,,ix'^y + 04,02;* 1 

y = x + 60, 2Z/^ + bi,ixy + h2flX^ + hi^2xy^ + b2,ix^y + &3,oa;'^ + ^2,22;^?;^ + b^^x^y + 64,0a;* J 

under the condition ( 12.61 ). That means finding all values of the 15 parameters for which the equation 
g'{x) + f{x) g{x) = 1 is satisfied where / and g are defined by (12.5b (see Corollary [OJ. 

When the standard reduction is possible, that means 04,0 = 03.0 — ^2,0 — 61.1 = 62,1 = 
63,1 =0 all the 6 isochronous centers with zero Urabe function were described in Theorem 3.1 
of ID. Otherwise when the Choudhury-Guha reduction needs to be applied the problem becomes 
substantially more complicated. 

Taking in account the great complexity of the problem we did not succeed in solving it com- 
pletely. Nevertheless, during our investigations we obtained 25 new isochronous centers for the sys- 
tem dl.llb . two of them of extreme complexity, called Monsters. We are convinced that our list is 
not exhaustive. 

The procedure to obtain the isochronous centers listed below consists in solving by Grobner 
method the system ( 12.61 ) simultaneously with the set of equations on parameters which corresponds 
the equation g'{x) + f{x) g{x) = 1. First, one applies the Solve routine of Maple (based on Grobner 
basis technic) which spHts the variety of solutions into 37 subvarieties. The cases ( I4.6b -( l4.23b were 
obtained by detailed inspection of some of them. The remaining 7 isochronous centers ( 14. 24b -( 14. 3 Ob 
were obtained by restricting ourselves to 62.2 = 03,1 =0 and by application of the standard Grobner 
basis technique. 
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We verified also that all above isochronous centers are not time-reversible, except perhaps the 
two Monsters ( 14.231 ) and ( 14.301 ). 

Theorem 4.2. The following quartic systems have an isochronous center at the origin O with zero 
Urate function. 



X = -y + ba^2xy + x^ - bo,2X^ 

y = X + bo^2y^ -2xy + 2x^ - bo.2x'^ 



X = -y + xy - as.ox^ + as.ox^ 

y ^x + y^ + 2 as^Qxy + 2 a^^'^x^ - a^^o^x* 



(4.6) 



(4.7) 



X = —y + xy — a^flx'^ + a^^x^ 

2 2 2/^1 2\ 3 2 4^' ^^'^^ 

y^x + 3y +2a3fixy-x + 4a3,oa: y + I - + 2a3,o jx + 03,0 a; 



X = —y + xy — a^fix'^ + a^^x^ 

3x'^ / 1 \ 

y^x + 4y^ + 2 a^^^xy - ^- + 6 a:iflx'^y + (l + 2 a3,o^) x"^ +[-- + 2 a^/ j x^ 

' (4.9) 



bo,2xy , 2 ^0,22; 



X = —y H h X 



3 3 

60,22^^ 4 6o,2a;^2/ , fbo^ , ^^ 3 , ^0,22;'* 
27 



4 ( , (4.10) 

y^x + bo.2y - 2 xy h — — + 2 x' 



, bQ^2xy 2 ^0.22;^ 

X = — y H ha: — 



, , 2 o 3 ^0.22:^ 3bo^2X^y (boa , ^\ 3 , f 1^3, ^0,2 \ 4 



(4.11) 



45 2 2 45 4 

x^-y-—xy + x +Y^ 

225 2 19 3 ., 3 ^' ^^'^^^ 

y = X — 2xy xy H a; + 45 x y 



x^ -y + 



J, 2 2 7 2 4 

64,0 a; y 2 04,0 a; 



2 2 ^^ (4 13) 

I J, 2 o O4 Qa; ,,22 oj, 2|o3 j,23ij, 4 

y = X + 64,02; - 2 a;y h 04,0 xy - 2 64,0a; y + 2x - 64,0 a; y + 64,0a; 
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y + f-2 + 2 VT9) x^y + x^ - (-2 + 2 VTd) x 



2 o aix 

y:^x±aiy -2xy^ ——- 



(4.14) 



Oi2xy T 2 aix y + a^x + a^x y ± 4 aix 



where 



a 



1 = \/-106 + 34a/T9, aa = -10 + 10 \/T9, ag = 16 - 16 \/l9, a^ = -13 + 3 719. 



22 2 s 4 2 

a; = —y + aiixy -\ ai.i x y + a; — .t ai^i x oi^i 



«i,iy^ o , 3ai,ix2 15 22,0 2,03^-- 23 

y = x 2xj/H h ^oi.i xy -|-3ai,ia; y + 2x - ^ ai,i x y 



15 9 q Sfli^ia;'* 
(4.15) 



x = -y^— a^xy + agx^y + x^ ± — a^x^ - a^x^ 



y = x± 



35 



2 xy =F — Q!5X^ + 5 a^xy"^ =F — a^x^y + a^x^ — 8 ag 2;^j/ ± — asx 



70 



35 



where 



as = V-77798+1162V469T, ae = -^ + ^^^ , ay 



35 



'^ + i^«l■ 



(4.16) 



Sas.ox-* 3 03,02;^ 
X = -y + xy h 03^02; — 



2 , 3a3,oxy 3 , Sag^ox^y 
y = X + 3 y + X + 



3 + g«3,o IX 



2^4 



19 2\ 3 303,02; y 3 03,0 X 



(4.17) 



v 2 2 2 V 2 3 V 2 4 
x = -y + xy±— -X =F^— X ±-^2; 
2 3 6 



y = X + 6 y =F v2xy x =p — v2x y - 



2 2 



13 3 , 3V2 , 

— X ± X y X 

3 2^3 



,3„, 4 4 



(4.18) 



X = — y T 2 a/2: 



y =F 2 V 2xy + x* 



±2v/2a 



14 



y = X T 6 V2y2 - 2 xy ± 2 V2x^ ± 8 V2x'^y + — x^ =F 2 V2x'' 

o 



(4.19) 



X = — y + fli^ixy + ai^i^x^y + x^ — ai^ix^ 

2 4 

— ai^i X 

y — x + 3ai,iy — 2xy — ai.ix +2ai.i xy 

— 4ai.ix y + 2x— 2ai,i x y + ai_ix 



(4.20) 
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2^4 



X ~ —y + ai^ixy + 3ai^i x y + x — ai_ia; — 3 Oi^i x 

y = X + Aai^iy — 2xy — — ai^x +6ai.i xy — Gai^ix y + 2x — 6aii x y + 2ai,ix 



(4.21) 



X = -y + ai^ixy + ai^i 



2 '^7 7 2 \ 2 



2 "^ I 2*^7 72\4 

+ 0; -aija; - ( ai,i -- &o, 201,1 + ^ °o.2 ) a; 



y = a; + 60,22/^ - 2 xy + ( - ai j - - 60,2 j a;^ + (2 oi^i^ - 3 60,201,1 + 60,2^) xy'^ 
2 f ai,i^ - - 60,201,1 + X 60,2^ ) ^ 2 (2 ai4^ - 3 60,201,1 + 60,2^) j x^y 



} ■ (4.22) 



+ (-2 ai,i + 60,2) x^ + (2 ai,i - 2 60,2) x'^y + 2x^ 



x = -y + Txy + Mx^y + .t^ - Tx'^ - Mx 
P 



y 



Py^ -2xy x'^ + 5 Mxy^ - 2 Px^y + 5x^-8 Mx^y + 4 Px'^ 



(4.23) 



where T = _RooiO/ (8436480 .^** - 151490752 .Z'' + 799019220 .Z'^ ~ 927412425 .Z'^ - 1728684180) 
has two real values and M = TJooiO/ (54647240223827376 .^^+ 

(561508456956838398 + 111131520938630016 T^ - 25008718903435904 T* + 1701065448975360 r«) .^^+ 
(17547249759911040 T^ - 5911680053099923488 + 2296193094798110772 T^ - 236397331187044016 T^) .Z''*+ 
(-47291962375532160 T** - 32577228112841141760 T^ + 1824502903342225392 T" - 38028219495361710336) _^^+ 
(54944501282517005964 - 141937291355478113601 T^ + 26240570003047890492 T^ - 3151542249299751840 T**) _^+ 
464619529558904831820+500960459979496602351 r2+20184288412145753196r'*-11100198119724668448T'^) 
has three real values. 

1277741426613483475261465072750940888801131320332201 \,T'irp , 387891696411877338314673950037666400855942884617917 f./j'^rp 

600834942501727375359564182536082155306896 -'"'"8718020422733426104722414152817729353615120239478544 -"^ -* 

1388592164504587266774323327561430923 a j^ 742426510189441463550540802769344501623332700250607 »J-37n3 

■391797844837578270259822642228874404 "'^ ^ 58846637853450626206876295531519673136902061616480172 -"^ -^ 



p _ 1277741426613483475261465072750940888801131320332201 \/l'irp , 38789169641187733831467395003766640085594288461791 
-" 78462183804600834942501727375359564182536082155306896 """8718020422733426104722414152817729353615120239478544 *- ' 
279706937387076388592164504587266774323327561430923 a j^ 742426510189441463550540802769344501623332700250607 »J-37n3 
242167233964817391797844837578270259822642228874404 "'^ ^ 58846637853450626206876295531519673136902061616480172 -"^ -^ 

1946020041492661760957343501655785409975201353232 M^T' I 525193640940998707118304158017605885033259029441052 )\/f3Tn 
14711659463362656551719073882879918284225515404120043 """220674891950439848275786108243198774263382731061800645 

103442863691365219385975139449969928492687792379916 a^2fTi5 164131380112569086879618812047151535419534003325851 nf2^3 
24519432438937760919531789804799863807042525673533405 6538515317050069578541810614613297015211340179608908 



lU4bU^UU414U^bblYbU9aY34daUlbaaY«a4UUUYa^U 13133^3^ nfo^' I b^blU3b4UU4UUUg<'UYllgd04ia«Ul<^b0b«gaU33^aUU^9441ua^ ?\, 
14711659463362656551719073882879918284225515404120043 """220674891950439848275786108243198774263382731061800645 

I 103442863691365219385975139449969928492687792379916 A/T^T^ 164131380112569086879618812047151535419534003325851 /iJ-2 
"""24519432438937760919531789804799863807042525673533405 6538515317050069578541810614613297015211340179608908 

367033244357164535298438083037260851819189099856 if2fT-,7 , 14771752963262353720563034395097445259447260691808 M^T' 
1634628829262517394635452653653324253802835044902227 """220674891950439848275786108243198774263382731061800645 

54666403014280865159557614303552794593967705728 »^4rp7 5151178974136622357147584757325383291000410195430 jj-- 



-.-UCi44ZgbCibbJlCibDZlU,JgOU JDlJl^44^^0bJ^^54bJ^0g ^ ^^^,3 ^^^Ib ^JArpO 1041Cil,3gUl IZOO^jUgOa ^ ^b Igg 1 ZU4 ^ 1 DIP J041^D J4UU J J^pgO 1 Ty^ArpO 

24519432438937760919531789804799863807042525673533405 6538515317050069578541810614613297015211340179608908 

367033244357164535298438083037260851819189099856 if2fT-,7 , 14771752963262353720563034395097445259447260691808 A/f^'TS 
1634628829262517394635452653653324253802835044902227 """220674891950439848275786108243198774263382731061800645 

54666403014280865159557614303552794593967705728 ;if4fTi7 5151178974136622357147584757325383291000410195430 M*T^ 
1711659463362656551719073882879918284225515404120043 14711659463362656551719073882879918284225515404120043 

4690601171400748048835548400017315870720 7n7 _ 121307988379654758804095121568312112060864 rj^5 

18911827251898810049102811291456652930124403 
I 36899089657974060583725199394228563447167805: 



' 14711659463362656551719073882879918284225515404120043 -- - 14Yllb5y4b33b2b5e55iri9Ur3K»:iKr99H 
4690601171400748048835548400017315870720 n^7 _ 121307988379654758804095121568312112060864 n^5 
"""6303942417299603349700937097152217643374801 18911827251898810049102811291456652930124403 

759721576540272607925488426699587086173816 ^"3 , 3689908965797406058372519939422856344716780525509 rp Jlfi 
6303942417299603349700937097152217643374801 """9807772975575104367812715921919945522817010269413362 
I 10767883711615521687536020314030171163247140 rp 82822172756429241023335015771633450343406027311056 ji/frpS 
""" 2101314139099867783233645699050739214458267 544876276420839131545150884551108084600945014967409 

152574332389107076817431693266653145437149456804321 \frp3 , 1501312340420936913937036728933692151818809426880 Ji/f^^ 
181625425473613043848383628183702694866981671655803 "*" 181625425473613043848383628183702694866981671655803 



M^ 



I 10767883711615521687536020314030171163247140 ^^ 82822172756429241023335015771633450343406027311056 
"•" 2101314139099867783233645699050739214458267 544876276420839131545150884551108084600945014967409 

I 152574332389107076817431693266653145437149456804321 \frp3 , 1501312340420936913937036728933692151818 
""" 181625425473613043848383628183702694866981671655803 """ 18162542547361304384838362818370269486698 

o _ 425588095896001129751500929401062088207356 /w-4 _ 5767486187745253514611173391155686944173385 
•-^ 510619335801267871325775904869329629113358881 1021238671602535742651551809738659258226717762 
I 2663265122906318241956827533449698379402528 _ 59835284352308243659547352890456360623975739 i#2 
""" 700438046366622594411215233016913071486089 226941927022785720589233735497479835161492836 
I 4276254433644147165491110877451438697226467 i^ , 20238619233578301504790424411472589676640 T-.6 
"""4202628278199735566467291398101478428916534 """ 700438046366622594411215233016913071486089 

193249058167562967816719266780463563218080 ^"6 »^ , 49806639391130680810149977028256722235840 T^6 a j-3 
18911827251898810049102811291456652930124403 """ 1531858007403803613977327714607988887340076643 

6528982704325582783377044452994315443840 fT-.6 a j2 , 5956207676950329133758464677648797191680 T^ A//"* 
170206445267089290441925301623109876371119627 """1531858007403803613977327714607988887340076643 
_ 687325278004245713061424259642842230089632 rp4 jw-2 _ 2050583245437666120312007255881326007529568 rp4 j./r3 

510619335801267871325775904869329629113358881 4595574022211410841931983143823966662020229929 

I 9808758961750027034269880327684876902994056 rpA t^/t 242134670257750264203544183229519465935616 T^4 /w-4 

""56735481755696430147308433874369958790373209 4595574022211410841931983143823966662020229929 
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_ 396343814205392842122901681054670066164296 rp4 , 2035287795094287693870135081243409761340961 rp2 » <-3 

700438046366622594411215233016913071486089 ^ "•" 1531858007403803613977327714607988887340076643 ^ ^" 
I 123744077281693368649685336016663155804936 T^ \f^ l 2168606272638745909274782927486819511187308 T^a iij-2 
"•" 1531858007403803613977327714607988887340076643 "*" 56735481755696430147308433874369958790373209 



14909180233466477515674693647959213932639919 nn2 33628313797235830289730557740669826423944793 rp2 



4202628278199735566467291398101478428916534 



37823654503797620098205622582913305860248806 



rM. 



X = —y + ai^ixy + x^ — ai,ix^ 



y = x + Aai^iy ~2xy- -ai,ix 



6 ai,ix^y + (2 + ai,/) .t^ + f 2 ai,i - - ax,i^ j 



3^x4 



(4.24) 



X — —y + a\,\xy + a; — a\^\x 

y = X + 3 a\,\y ^ 2 xy — ai,ix'^ — 4 ai^ix'^y + ( — ai.i"^ + 2 ) .t^ + aiix'' 



(4.25) 



i = -y + 2 /3 xy + .T^ - 2 /3 x^ 

y = a; + 8/3 2/2 - 2x^-3 /3a;2- 12 ;3a;2y+ 14x^-2 /3 a;'' 



where /3 = ± a/3- 



X = — 2/ + axj/ + x axH — x 



y = x + 3ay — 2xy — ax — 3axyH — x+2 x y ax 



where a = ± \/2- 



(4.26) 



X = — y + a xy + X — 4/3 a x + 2/3 x 

, 2 2 26 , , , 4^ (4.27) 

i/ = x + 6ay — 2xy — 5/2ax — 9Q;xy+— x+6 x y — 8/3 a x 



(4.28) 



X = -y + aisxy + x^ + ( oi^i^ - - 60,201,1 + 2 ^0,2^ ] x^V - ai^ix^ + ( -aij^ + - 60,201,1 - 2 ^0,2^ ] a;"* 

y = X + &o,2y^ - 2 xy + f - oi^ - - 60,2 j x^ + (2 ai,i^ - 3 60,201,1 + 60,2^) xy"^ 

+ (2 fli.i - 2 60,2) x^y + 2x^ + (-2 a^^ + 3 60.201,1 - 60,2^) x^y + (-2 ai.i + 60.2) x^ 

(4.29) 
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X = —y + axy + M X y + X — ax — M x 

S 2 -^3.0 3 -84,0 4 

-X y ^o ^-i 



(4.30) 



3{M + 9)af3 " 2{M + 9)f3 12(M + 9)a/3 

where (M + 9)af3^0 and a is the real root of the polynomial P{Z) : 

P{Z) = (144 M + 1296) .Z*^ + (7236 M - 24 M^ + 996 M^ - 15552) _Z^ 

+ (12960 M2 + 2736 M^ - 202 M^ - 46656 M) .Z^ 

-8 M6 + 664 M^ - 109 M^ - 15552 M^ - 233280 M + 10368 M^ 

and 

13 = 1581568416-224710848 M3+452213280 M2-20366712 M4+284739 M6+3438793 M'^+ 

32 M^ - 892 M« - 15958 M^ + 2301021432 a^ + 655618860 a^ + 4745755008 M + 96 a*M'^ - 

40132 a^M^ + 1000 M^a^ _ 3^20 a^M^ + 128736 a^M* + 530942 a^M^ - 178814952 M ^a^ + 

4111938 M4a2+1440585648 ^2^-68912208 M3a2-36318780 Ma^+216852 M^a^ -22879908 M^a* 



S = - 1309805864856 M3+1304677734768 M^- 109564995222 M4+3102777711 M6+32534070309 M^+ 

288136 M9-31013013 M8-281741535 M^+44302 Mio+341618777856 a2+1640171477952 a*+ 

1000 ikfiOa^ + 9g Q,4^9 _ ;Lg8o a^M^ + 17574286201920 M+ 1910484 a^M^ - 49132 a^M^ + 

4176 ^4^^-4906584 a^M 6-286521084 a2M6-30232362 M^a2^34671168 a4M5+1921156704 a^M^+ 

7550896788 a2M5-708776248212 M2a2+32149733976 M4a2^4446970742160 a2M-333064034298 M^a^- 

226136884752 Ma* + 4800914064 M^a* - 99633193992 M^a* - 796 M" + 32 M^^, 



Bo,2 = 854046944640 0^+4378535313120 a4+18416851416 M3a''-263459876832 M'^a*- 
711485573796 Ma'' + 11974507973856 a^M - 884184180804 M^a^ + 104273686746 M^a^ _ 
2161040187528 M2a2_^20100360978 a2M5+5321931120 a4M4+30597468 a4M5-16021728 a^M^- 
1086777654 a^M 6-80729938 M V+93312 a''M^-179920 a2M9+6307872 a2M8-5568 a'^M^+ 
384 a4M9+4000 MiOa2-3677920975152 M ^-249464899332 M4+7497589581 M6+95584892967 M^+ 
140776 Mio+1979888 M9-79612855 M«-920925749 M^-2800 M"+128 Mi2+47893803042240 M+ 
2753000428896 M^ 

^3,0 = -56936462976-139773234528 a2-204261985872a4+14956272M3a4_^;L;L44gg48204M2a4+ 
13633897968 Ma4_4894g778g7^2 a2M+35663387130 M3a2-1670720580 M4a2^51974431164 MW- 
718427376 a^^^ - 165970176 a^M^ - 5669352 a^M^ + 260880 a^M^ + 8347234 a^M^ + 
2613748 M V - 3696 a^M^ + 3000 a^M^ - 119572 a^M^ + 288 a*M^ + 117311420664 M^ + 
14135146866 M* - 342595189 M^ - 2512013814 M^ - 1940 M^° + 36898 M^ + 2916172 M« + 
15602085 M^ + 96 M" - 1772624526048 M - 200707918416 M^, 



B4.0 = -512428166784 a^ + 1155128820192 a* + 76187990232 M V - 37449803088 M^a* - 
1560658377780 Ma* + 4471598426688 a^M - 170217184320 M^a^ + 262466972514 M*a^ - 
3995129998656 M2a2+4574201922 a2M5+3865776624 a4M4-781052628 a4M5-41458608 a^M^- 
4386392646 a^AfS - 21588886 MV + 1071072 a*M^ - 668472 a^^^ + 17045880 a^M*^ - 
15264 a4M8+1728 a''M9+18000 M iOa2_327957327i872 M3+482166146232 M4-8025777621 M'^+ 
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136445218065 M5+323820 Mi°+15941920 M9+19475327 MS-2392822547 M^-9336 M^^~ 
576 Mi2 + 25093708683840 M - 8570037406464 M^. 



5. Explicit Linearization 

5.1. Linearization formulas 

Let us consider the Lienard type system ( 11.4b 

x = y 1 

y = -g{x)- f{x)y'^] 

with a center at the origin (0, 0) where / and g are real analytic in a neighborhood of zero. 

It is known by Sabatini formula ( 11.6b that the first integral associated to the system ( 11.4b can 
be written 



I(x,x)^ f g{s)e^^<-'U.s + liie^'-'^^y' 



2- 

where F{x) = J^ f{s)ds. 

Following 121 (see also |fT9i ). let us perform the following change of variables 

/ . \ . F(x) 

p[x, x) = xe ^ ' 



(5.1) 



(5.2) 



As -——^ — - — — e^-'^(^) < and p(0, 0) = q(0) = then this is a change of variables 
o{x,xj 

preserving the origin and well defined around it. Moreover, q'{x) = e^^^'> > and thus the function 

X I — > q{x) is strictly increasing. In the {p, q) coordinates the first integral (15.11) becomes 

I{x,x)^H{p,q)^^p^ + Uiq), (5.3) 

where U is some uniquely defined real analytic function, U{0) = 0. Now it is easy to see that the 
system ( 11.4b in {p, q) coordinates can be written as 



dH 

q^ ^T ^P 

op 

p^-^^^±U 
dq dq 

that is as a Hamiltonian system corresponding to the Hamiltonian ( 15.31 ). 

The main result of this Section is 



(5.4) 



Theorem 5.1. Let us consider the Lienard type system ( 11.41) with real analytic functions f and g 
such that X g{x) > for x ^ Q. Then the origin O is an isochronous center with Urate function 

2 

ft, = if and only ifU{q) = ^. 

Proof Taking into account our assumptions, from Theorem ll.il one knows that O is a center. Now, 
from Corollary 11.31 one knows that O is an isochronous center with Urabe function h — Q if and 
only if (7'(x) + g{x)f{x) = 1 or equivalently f;'(x)e^*^^^ + g{x)f{x)e^^'^^ — e-'^^^h The last equal- 
ity is nothing else (^g{x)e^'^^^^ = (/„^ e^'-''^ds) . As (7(0)e^'-'^' — 0, when integrating one ob- 
tains g{x)e^''^^ = J^ e-'^^^^ds or equivalently U'{q) ~ q because {^){q{x)) — g{x) e^'^^K Since 
C/(0) = Gone has C/(q) = ig^. D 
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Consequently when the Urabe function h identically vanishes , the system of coordinate (p, q) 
defined by ( I5.2l i is the linearizing system of coordinates for system ( ll.4l i. Indeed, 



q=p 
p = -q 



(5.5) 



It is interesting to compare the theorem ISTI with Chalykh-Veselov theorem that we formulate 
only for potential U without pole at 0. 

Theorem 5.2 (|4|, Theorem 1). Let us consider the Hamiltonian system with the Hamiltonian 
H{P, Q) = iP'^ + U{Q) where U is a rational function without pole at 0. Then O is an isochronous 

rifi rifi 

center for the associated Hamiltonian system Q = — — , P = — — — if and only if up to a shift 

oP oQ 

Q ^)- Q + a and adding a constant, U{Q) — k Q^ for some fc e M — {0}. 

5.2. Examples 

Now applying formula (15.2b we provide 5 examples of linearization of isochronous centers with 
zero Urabe function. The reduction to Lienard type system (ll.4l i is always obtained by standard or 
Choudhury-Guha reduction. To compute variables {p, q) see ( (I5.2l i) we use Maple, and the identity 
( 15.5b was verified in all cases. Our choice is somewhere random, because all reported examples with 
zero Urabe function are good for such purpose. 



5.2.1. Cubic examples. 

1 . Consider the case 6 of Theorem 3 of lfT3l . that is the system 

, , ( 2 3aii6o,2 , ^0,2 \ 2 

X = -y + ai^ixy + \ais ^ \ ^j^V 

y ^ x+ i — y- + -i;i j x'^ +bo,2y^ + (2 Oi^i^ - 3ai_i6o,2 + ^0,2^) xy'^ 
In this case the functions / and g are 



(5.6) 



/(^) 



^0.; 



(2ai,i2 - 3 ai,i6o,2 + 60,2') x + ai,i + 2 {a,,i^ - ^^^i^ 



bo,: 



-1 + aiax + (ai,i' ~ ^""2'°'" + ^) ^^ 



9{x) 



1 — ai,ix 



2 3 ai 160.2 , ^0,2 \ 2 
ai.i 1 ^r- I X 



2 2 

for which we obtain the following linearizing change of coordinates 



^0.2 , ai,i \ 2 



(A(x)-A(O)) 



q[x) 

p{x,y) 



2 (01,1-60,2)^ 

(2 + ai^ix - 60,2a;) xe V-""i>i'+''"i-i''o,2-2t.o,22 

(—2 + 2 ai^ix + 2 oi.i^x^ — 3 x^ai_i6o,2 + x'^bo,2 ) 
-2yq{x) 



where 



(2 + ai_ix — 60,2a;) X 

2 oi.i^x — 3 a; oi 160,2 + a;6o,2^ + ai, 



>, 



(5.7) 



(5.8) 



A{x) = arctan 



-5 ai^i^ + 6 ai,i6o,2 — 2 60,2 



(5.9) 
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2. Consider system (I4.3l l from Theorem l4.1l In this case the functions / and g are 



/(^) 



^2,0 



1 + 262,02; 

g{x) = a; (2 foa.o^aj^ + 3 62,0a; + l) 
for which we obtain the following linearizing change of coordinates 

.X (62,02; + 1) "1 



(5.10) 



q{x) 



p{x,y) 



(1 + 262,00;)' 

— y + x'^ 
(1 + 262,02;)' J 



(5.11) 



5.2.2. Quartic examples. 

1 . As a quartic example we consider the system III of Theorem 3.1 of f3l . We choose the follow- 
ing restrictions on the parameters 021 = —3, 60,2 = —4, 620 = 5 to obtain simple, presentable 
expressions for linearizing variables. 



X = —y — 3xy — 3x y — x y 

y = X + — X — 4y — 4xy — 2x y 

In this case the functions / and g are 

7 + 10a; + 5a;2 



(5.12) 



/(^) 



l + 3a; + 3a:2 + a;3 



1 



(5.13) 



g{x) = - {l + 3x + 3x'^ + x^)x{2 + x) 
for which we obtain the following linearizing change of coordinates 



1 _ , _,;,>^^ ^^ ^ ^^^_2 



q(x) ^ ^x{2 + x)e (1+-)^ (1 + 2;) 



(5.14) 



Pix, y) 



-ye 



{i + xY 



2. We consider the system ( 14.231 ) of Theorem l4.2l 
In this case the functions / and g are 

., , 4 ai,i^a; - 6x60,201,1 + 2x60,2' + ail +60,2 

j{x) — ■ ■ — 2 

1 — ai,i2; — ai,i^a;^ + 3/2 a;^6o, 2^1,1 — 1/2 2;^6o,2 

g{x) — —1/4 X (ai,ia; — a;6o,2 + 2) (—2 + 2 ai,ia; + 2 ai,i^a:;^ — 3 a;^6o, 2^1,1 + 2;^6o,2 ) _ 



for which we obtain the following linearizing change of coordinates 

X (fli.ix — a:;6o,2 + 2) 5 (x) 



q{x) = - 
p{x,y) = 2 



where 



2(<il,l~i'0,2) arct 



—2 + 2 ai,ia; + 2 ai^i^x^ — 3 2;26o, 201,1 + x^bo,2 

jx' 'y)S jx) 

—2 + 2 ai,ia; + 2 ai,i^a;^ — 3 2;^6o, 201,1 + 2:^60,2 , 



2aii 3J-3 a:bo, 2'^!, 1+^^0,2 +^1,1 



> , 



(5.15) 



(5.16) 



-\/-S"l,l^+8i>o,2"l,l-2bo,2^ 



■y/-5°l,l^+6bo,2fl,l-2bo,2^ 



S{x) 



-^-Soi. 1^+6 to. 2 "1,1 -2 bo. 2^ 
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5.2.3. A rational example. Let us consider the system 

yx 



X = -y 
y = x + 



l + x , 
2 > ■ (5.17) 



y 



In this case the functions / and g are 



1 + x 

2 + x 



f>, (5.18) 



gi^) 



l + x 
for which we obtain the following Unearizing change of coordinates 

q(x) = xe^ 1 

f ^ -n , ^\■ (5.19) 

Pix,y) =2/e="(l + a;)J 

5.3. Comments 

It is really astonishing that in all above cases the linearizing variables {p, q) are always expressed in 
"finite terms". This follows from the fact that if ^'(a;) + f{x) g{x) = 1 then / — ^—^. Moreover, 
as g{0) = and ^'(O) = 1 the singularity of / at zero is spurious. In all examples considered in 
this and related papers lfT2l[T3l [3]| / and g always are rational functions. Then F{x) — J^ f{s) ds is 
expressed in "finite terms" and thus also p{x, x). The problem is slightly more delicate in what con- 
cerns q{x). But / e^(^) ds = j e^ ^^^^ '^^ ds = J e^ ^^^^ "^^ ds ^ J -^e^ d^ '^' ds ^ e^ '^ ''" + 
const, g being a rational function, J -4jr is obtained in "finite terms" and thus also J e^^*-' ds. 
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